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Abstract 

We have studied the effects of momentum dependent interactions on the single-particle 
properties of hot asymmetric nuclear matter. In particular, the single-particle potential of 
protons and neutrons as well as the symmetry potential have been studied within a self- 
consistent model using a momentum dependent effective interaction. In addition, the isospin 
splitting of the effective mass has been derived from the above model. In each case tempera- 
ture effects have been included and analyzed. The role of the specific parametrization of the 
effective interaction used in the present work has been investigated. It has been concluded 
that the behavior of the symmetry potential depends strongly on the parametrization of 
the interaction part of the energy density and the momentum dependence of the regulator 
function. The effects of the parametrization have been found to be less pronounced on the 
isospin mass splitting. 

PACS number(s): 21.65.+f, 21.30.Fe, 24.10.Pa, 26.60. +c 

1 Introduction 

One of the most interesting problems in nuclear physics is the isovector dependence on nuclear 
force, which can be found in nuclear symmetry energy, the isovector optical potential and neutron- 
proton effective mass splitting. The isovector feature of nuclear forces is crucial in order to gain 
a good understanding of neutron stars and exotic nuclear collisions produced at radioactive beam 
facilities and to describe the structure of exotic nuclei. This interest in the isospin dependence 
of nuclear forces is of recent date because data for neutron-rich nuclei were rather scarce in the 
past. The forthcoming new generation of radioactive beam facilities, such as the future GSI 
facility FAIR, the Rare Isotope Accelerator planned in the USA and the SPIRAL2 at GANIL, will 
produce huge amounts of new data for neutron-rich nuclei. Up to now, the isovector dependence 
of nuclear force has been investigated in heavy ion collisions experiments. The advantage of such 
kinds of reactions is that they allow testing of nuclear forces at supranormal densities since in 
intermediate energy compressions of two to three times nuclear saturation density is reached. 
Nevertheless, the asymmetry of the colliding systems is moderate and therefore the isospin effects 
on the corresponding observables are moderate as well. 

From a theoretical point of view, the predictions for the isospin dependence of nuclear interac- 
tion are very different. In general, both microscopic and effective interactions have been extensively 



1 



used in order to gain knowledge about the nuclear matter properties in conditions far from equi- 
librium (hot nuclear matter, high density behavior and so on). It is appropriate therefore, in every 
case to incorporate temperature and momentum dependence in order to have a richer interaction 
and as a consequence be able to produce a more thorough description of nuclear matter properties.. 
Specifically, the determination of the nuclear symmetry energy (NSE) based on microscopic and/or 
phenomenological approaches is of great interest to nuclear physics and nuclear astrophysics. For 
instance, it is important for the study of the structure and reactions of neutron-rich nuclei, the 
Type II supernova explosions, neutron-star mergers and the stability of neutron stars. So far, 
the main part of the calculations concerning the density dependence of the SE is related to cold 
nuclear matter (T = 0). However, recently there has been an increasing interest in the study of 
the properties of asymmetric nuclear matter, including NSE, and the properties of neutron stars 
at finite temperature dDl [HI [IS [151 ESI dZl [IBl HSl 1201 121] 

In our previous work [21] we studied the effects of finite temperature on NSE and we also 
found the appropriate relations describing that effect. We especially focused on the interaction 
part of the NSE, which so far has received little theoretical attention concerning its dependence 
on temperature. We applied a momentum dependent effective interaction model. This model, 
afterwards called BGBD (Bombaci, Gale, Bertsch, Das Gupta), has been introduced by Gale et al. 
[221 [231 [2H [22] to examine the influence of momentum-dependent interactions on the collective 
flow of heavy ion collisions. Over the years the model has been extensively applied to study not 
only the heavy ion collisions, but also the properties of nuclear matter by proper modification 
[H [2H1 [23 12S 129]. A review analysis of the present model is presented in Refs. [U |2"3] . 

In the present work, we have studied the momentum and temperature dependence of the mean 
field properties of asymmetric nuclear matter |^|Sll^l^Ell^l^EZll3Ell3SlllDlinill2lll31 
SH S3, SSI S71 SHI S9l [501 EH E2]- Efforts up to now have been devoted mostly to studying the 
properties of cold nuclear matter. In contrast, the motivation for the present work is to study 
the properties of hot asymmetric nuclear matter, especially the temperature dependence of the 
nuclear symmetry energy and the single particle properties of nuclear matter. 

We have employed a model with the characteristic property that the interacting part of the 
energy density is momentum dependent. Thus, the single particle potentials of protons, neutrons, 
as well as the symmetry potential, are also momentum dependent. The isovector part of the 
optical potential, i.e. the symmetry potential, describes the difference between the neutron and 
proton single particle potentials in neutron rich matter. The symmetry potential is one of the 
basic inputs to the transport models for the collisions of radioactive nuclei [30]. In addition, due 
to the momentum dependence, the temperature is expected to affect not only the kinetic part, 
but also the interacting part of the energy density. This is important in the sense that the density 
dependence of the nuclear symmetry energy, influenced by temperature, has a powerful effect on 
the values of the proton fraction and consequently the composition of hot /3-stable nuclear matter, 
with extensive applications in heavy-ion collisions and Nuclear Astrophysics. 

The effective mass (EM) is one of the most fundamental properties characterizing the propa- 
gation of a nucleon in a nuclear medium. Knowledge about nucleon EM in neutron rich matter is 
crucial to fully understand several properties of neutron stars. EM is determined by the derivative 
of the single particle potential with respect to the momenta k for k = hp. Thus, the trend of 
the EM is directly connected to the momentum dependence of the corresponding single particle 
potential. Furthermore, the isospin splitting of the effective mass i.e. the difference between the 
neutron and proton effective masses is derived from the above model. The present study is a 
contribution to the theoretical study of the neutron-proton effective mass splitting, a problem 
which is still highly controversial within different approaches and/or using different nuclear effec- 
tive interactions [5]. The above analysis indicates the necessity to apply a momentum dependent 
interacting model to study the single particle properties of hot nuclear matter. 
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This work is a continuation of recent papers [23 E2] , where the authors have employed a phe- 
nomenological non-relativistic effective interaction first introduced in Refs. [H[26], with suitable 
modification concerning the isovector part of the interaction. Here, we have applied a more gener- 
alized expression of the interaction part of the energy density, with a richer parametrization and 
additional parameters introduced to maintain causality pQ . We have concentrated on a systematic 
study of the effect of the parametrization of the effective interaction on the mean field properties 
of hot asymmetric nuclear matter. 

The plan of the paper is as follows. In Sec. II the model and the related formulae are discussed 
and analyzed. Results are reported and discussed in Sec. Ill, while Sec. IV contains a summary. 

2 The model 



The schematic potential model used in the present work, is designed to reproduce the results of 
microscopic calculations of both nuclear and neutron-rich matter at zero temperature and can 
be extended to finite temperature [H [2S] • The energy density of the asymmetric nuclear matter 
(ANM) is given by the relation 



e(n n , rip, T) = e% in (n n , T) + e p kin (n p , T) + V int {n n , n p , T), 



where n n (n p ) is the neutron (proton) density and the total baryon density is n = n n + n p . The 
contribution of the kinetic parts is 



d 3 k h 2 k 2 

4in( n n, T) + e p kin (n pj T) = 2 I (/«K, K T) + f p (n p , k, T)) 



(2tt) 3 2m 

where f T , (for t = n,p) is the Fermi-Dirac distribution function with the form 

'e T (n T , k,T) - /i r (n r ,T) N 1 



(2) 



f T (n T ,k,T) 



1 + exp 



T 



(3) 



This distribution is inserted into the following integral in order to evaluate the nucleon density 

TV, 

d 3 k 



1l T 



(27T) 



-f T (n T ,k,T). 



In Eq. ([3]), e T (n T , k,T) is the single particle energy (SPE) having the form 



e r (n r , k,T) 



2 1.3 



h k 
2m 



+ U T (n T , k,T). 



(4) 



(5) 



fi T (n T ,T) stands for the chemical potential of each species, while the single particle potential 
U T (n T , k,T) is obtained by the functional derivative of the interaction part of the energy density 
with respect to the distribution function f T . 

Including the effect of finite-range forces between nucleons to avoid acausal behavior at high 
densities, the potential contribution is parameterized as follows 



V in t(n n ,n p ,T) 
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where 

J>,I,T)= 2 J ^LgfaAjfT- (7) 

In Eq. ([6]), I = (n n — n p )/n and u = n/no, with uq denoting the equilibrium symmetric nuclear 
matter density no = 0.16 fm~ 3 . The parameters A, B, a, C\, C% and B', which appear in the 
description of symmetric nuclear matter and the additional parameters Xq, £3, Z\, and Z 2 used to 
determine the properties of asymmetric nuclear matter, are treated as parameters constrained by 
empirical knowledge [H [26] . 

The first two terms of the right-hand side of Eq. ([6]) arise from local contact nuclear interaction 
which led to power density contributions such as in the standard Skyrme equation of state. These 
are assumed to be temperature independent. The third term describes the effects of finite range 
interactions, according to the choice of the function g(k,Ai), and is the temperature dependent 
part of the interaction [26]. The function, g(k,Ai), may have the following forms 

1. Case 1: 



gi(k,Ai) = 




In this case we introduce two finite range terms: one corresponding to a long-range attraction 
and the other to a short-range repulsion. The finite range parameters are Ai = 1.5k F and 
A 2 = 3k F and k° F is the Fermi momentum at the saturation point n Q . The function g±(k, Aj) 
has been used extensively in previous papers (see Refs. [H EH] and references therein). 

2. Case 2: 



g 2 {k,Ai) = 




In this case the finite range interactions are approximated by effective local interactions 
by retaining only the quadratic momentum dependence. Therefore, the energy density in 
Eq. ([I]) takes the form of Skyrme's effective interactions. Actually, as we will show later, the 
two functions coincide, for low value of momenta k (k < 1 fm _1 ), but they exhibit different 
trends for high values of k. 

The energy density of asymmetric nuclear matter at density n and temperature T, in good 
approximation, is expressed as 

e(n, T, I) = e(n, T, I = 0) + e sym {n, T, I), (10) 

where 

e sym (n, T, I) = nI 2 Eif m (n, T) = nl 2 (jS^(n, T) + E™ m (n, T)) . (11) 

In Eq. (ITT1) the nuclear symmetry energy El^ m (n,T), is separated in two parts i.e. E^Jji, T) 
(kinetic) and E % ^ m {n, T) (interaction). 

From Eqs. (1101) and fill I) and setting I = 1, we obtain that the nuclear symmetry energy 

E ly m ( n i T ) is S iven b y 

E% m (n, T) — — (e(n, T,I — 1) — e(n, T, 1 = 0)). (12) 

Thus, from Eqs. (fl2"j) and (TjQ) and a suitable choice of the parameters xo, x%, Z\ and Z 2 , we can 
obtain different forms for the density dependence of the symmetry energy E^ m (n,T). It is well 
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known that the need to explore different forms for E^ m (n,T) stems from the uncertain behavior 
at high density pQ . The high-density behavior of symmetry energy is the least known property of 
dense matter [531 EH G2] , with different nuclear models giving contradictory predictions. Thus, 
in relativistic mean field (RMF) models, symmetry energy strongly increases with the density of 
nuclear matter [56], while in many realistic potential models of nuclear matter in the variational 
approach [31], the symmetry energy saturates and then bends over at higher densities. 

Recently, the density dependence of the symmetry energy in the equation of state of isospin 
asymmetric nuclear matter has been studied using isoscaling of the fragment yields and the an- 
tisymmetrized molecular dynamic calculation [57J. It was observed that the experimental data 
at low densities are consistent with the form of symmetry energy, E sym (u) ~ 31.6m ' 69 , in close 
agreement with those predicted by the results of variational many-body calculations. In Ref. [57J 
it was suggested also that the heavy ion studies favor a dependence of the form E sym {u) ~ 31.6m 7 , 
where 7 = 0.6 — 1.05. This constrains the form of the density dependence of the symmetry energy 
at higher densities, ruling out an extremely "stiff" and "soft" dependence |57j . 

Additionally, by using the isospin dependent Boltzmann-Uehling-Uhlenbeck transport model 
calculations, Chen et al. [29] also showed that a stiff density dependence of the symmetry energy 
parameterized as E sym (u) ~ 31.6m 105 clearly explains the isospin diffusion data [58j from NSCL- 
MSU (National Superconducting Cyclotron Laboratory at Michigan State University). 

In the present work, since we are interested mainly in the study of thermal effects on the nuclear 
symmetry energy, we choose a specific form for it, enabling us to accurately reproduce the results 
of many other theoretical studies [591 [60]. I n R- e f- [SH] the authors carried out a systematic analysis 
of the nuclear symmetry energy in the formalism of the relativistic Dirac-Brueckner-Hartree-Fock 
approach, using the Bonn one-boson-exchange potential. In a very recent work [60j, the authors 
applied a similar method to that in Ref. [59] for the microscopic predictions of the equation of 
state in asymmetric nuclear matter. In that case E sym (u) was obtained by employing the simple 
parametrization E sym (u) = Cu 1 with 7 = 0.7 — 1.0 and C ~ 32 MeV. The authors concluded that 
a value of 7 close to 0.8 gives a reasonable description of their predictions, although the use of 
different functions in different density regions may be best for an optimal fit [60]. The results of 
Refs. [59l [60] are well reproduced by parameterizing the nuclear symmetry energy according to 
the formula 

E% m (n,T = 0) = 13m^+ £F(m) . (13) 

Kinetic Interaction 

For the function F(u), which parametrizes the interaction part of the SE, we apply the following 
three different cases 

2m 2 

F 1 (m) = v ^I, F 2 (m)=m, F 3 (u) = — — . (14) 

1 + M 

The parameters x , x 3 , Zi and Z 2 are chosen so that Eq. f[T2|) . for T = 0, reproduces the results 
of Eq. f[T3|) for the three different forms of the function F(u). In addition, the parameters A, B, 
a, C\, C2 and B' are determined in order that E(n = n Q ) — mc 2 = —16 MeV, n = 0.16 fm -3 , 
and the incompressibility to be K = 120, 180, 240 MeV for each of the three cases. 



2.1 Single particle potentials 

The single particle energy e r , obtained by the functional derivative of the energy density (Eq. flTJ) 
with respect to the distribution function f T , is written as 

h 2 k 2 

e T {n,I,k,T) = + U T (n,I,k,T). (15) 
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The single particle energy e T consists of a kinetic part and an interaction one U T (n,I,k,T), 
which depends explicitly on density, momentum, isospin asymmetry and temperature as expected 
from an interaction term. 

The single particle potential U T (n, I, k, T) (protons or neutrons), obtained from the functional 
derivative of the interaction part of the energy density (Eq. (j^D) with respect to the distribution 
function f T , has the general form 



U T {n, I, k, T) = U?{n, I) + U?(n, I) + U^ D (n, I, k, T). 



(16) 



The first two terms are momentum independent, while the third one describes the momentum 
dependence of the single particle potential. The three terms have the following forms 



with 



U*(n,I) 



U?{n, I) = Aut ^A{^ + x )ul, 

U^(n, I)U? 2 (n, I) - U* 3 (n, I)U?*(n, I) 
[U T m (n,I)] 2 



(17) 
(18) 



and 



U T B1 (n,I) 
U T m (nJ) 
U T m (n,I) 
U T m (n,I) 



B{a + l)u a =f ^BQ + x 3 )u a I + ~B(1 - a)Q + x 3 )u a I 2 



1 + -B' 

3 

B '( 
— {a - 

n 
2 

3 Bn ° 



u 



0-1 



l)u^ T \-{\ + x,)u^I 
6 n 2 



^(3-.)(i 
3n 2 



x 3 u 



0-2 t-2 



U T MD (n,I,k,T) 



5 ^0 



u 



0+1 



(19) 



- (3C t - AZi) T T + (Q + 2Zi) J l T , 



+ «E 



i=l,2 



a± 



Ci — 8Zi 



I)g(k,Ai 



(20) 



The subscripts in the integrals are r ^ t'\ the upper signs stand for neutrons, while the lower 
ones for protons. An advantage of the present model is that for T = 0, the term U^ [D (n, I, k, T), 
as well as all the quantities can be derived in analytical forms. 

It is of interest to see that the single particle potentials are separated in two parts. The first 
one, U^(n,I) + U^(n,I) is a function only of the baryon density n and the isospin asymmetry 
parameter /. The second one, U^ 4D (n, I, k, T) has an additional dependence on T and k. Actually, 
U^ rD (n, I, k, T) is mainly responsible for the trend of the effective mass and also the effective mass 
splitting. Additionally, it is connected with the effect of the temperature on the interacting part 
of the energy density. It is also worthwhile to notice the direct correlation between the regulator 
function g(k, Aj) and U^ D {n, /, k, T). Thus, the choice of the function g(k, Aj) plays an important 
role for the single particle properties of hot nuclear matter, but one expects it to also be significant 
for the bulk properties of nuclear matter. 



2.2 Nuclear symmetry potential 

The nuclear symmetry potential (NSP) refers to the isovector part of the nucleon mean-field 
potential in isospin asymmetric nuclear matter, which in hot nuclear matter can also depend 
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on the temperature. Most of the studies concerning the NSP have been carried out for zero 
temperature, while the temperature dependence of the NSP so far has received little theoretical 
attention [8]. The NSP potential, describes the difference between the neutron and proton single 
particle potentials in neutron rich matter and has the form 



U, am (n,IXT) = U " {n ' I ' k ' T) - I U ' {n ''' k ' T) . (21) 



Various theoretical models have been applied to study the symmetry potential. Most of them 
predict a symmetry potential decreasing with increasing nucleon momentum. However, some 
nuclear models which employed effective interaction, predict an opposite behavior [8]. 

A systematic analysis of a large number of nucleon-nucleus scattering experiments and (p,n) 
charge-exchange reactions at beam energies up to 100 MeV has shown that the data can be very 
well described by the parametrization 

Usym^Ej^in) d (22) 

with a m 22 — 34 MeV and b = 0.1 — 0.2. Actually, the uncertainties for both parameters a and b 
are large. As pointed out in Ref. [15], the old analysis of Lane [61] is consistent with a decreasing 
trend of the potential as function of k, while a more recent analysis based on Dirac phenomenology 
[62] leads to the opposite conclusions. 

In order to clarify the effects of the momentum dependence on the NSP, it is easy to show, by 
applying Eqs. ffT6l) - fT20|) . that U sym (n,I,k,T), at zero temperature is 



U sym (n, I, k) = U^ p (n, I) + U^ p (n, k). (23) 



In Eq. (|23|) . 17^^ (n, I) is momentum independent, while the momentum dependent part U^® p (n, k) 
is written as 

U^ P (n, k) = - £ (C t - 8Z>) g(k, A,). (24) 

i=l,2 

From Eq. ( 124"1) it is obvious that the behavior of U sym (n, k) as a function of k is defined from the 
values of the parameters C{ and Zj, as well as from the function g(k,Ai). In fact, the regulator 
function g(k,Ai) specifies the trend of the slope. Thus, it is interesting to study the effect of 
the parametrization, including the function F(u), and the contribution of the choice of a specific 
function g(k,Ai) on the properties of the NSP, both for cold and hot nuclear matter. 



2.3 Effective mass 



The nucleon effective mass is one of the most important single-particle properties of nuclear 
matter. It characterizes the momentum dependence of the single-particle potential of a nucleon, 
and consequently the quasiparticle properties of a nucleon inside a strongly interacting medium 
such as the nuclear matter. Moreover, the effective mass describes to leading order the effects 
related to the nonlocality of the underlying nuclear effective interaction and the Pauli exchange 
effects in many-fermion systems |4"5| [371 [66] . 

The effective mass m*.(k), is determined by the momentum-dependent single nucleon potential 

via 



(n, J, k) 



1 + 



dU T (n, I, k, T) 



h 2 k 



dk 



(25) 



According to Eq. (1251) . m* (n, /, k) generally depends on the baryon density, the isospin asymmetry 
and the momentum of the nucleon. In the present model in* is independent of T. 



7 



An evaluation of m* mass at the Fermi momentum k = kp, employing Eq. (|25|) . yields the 
Landau effective mass. An advantage of the present model is that by applying Eq. ([161) . we get 
m* in analytical form i.e. 



m*(n,I) 



m. 



m 



u- 



h k 



F i=l,2 



dk 



I k=kp 



(26) 



Eq. fl26|) exhibits the dependence of the effective mass on the values of the parameters Cj, Zi and 
Aj as well as on the derivative of the regulator function g(k,Ai). In particular, for the two cases 
of the function g(k,Ai) (given in Eqs. (JSJ) , and Q), we get respectively for the effective masses 



m*(n, I) 



2wm T 



E 



i (a ± °*=*&i 



A 2 

i=l,2 



;..() 



l+ S ((i±/)«) 



2/3 



-1 



case 



ml 



(n,I) 



2um T . , 1 
^ A 2 

i=l,2 J H 



a 2 



case 



(27) 



(2f 



The meaning of Eqs. (1271) and ( 1281) is clear. Firstly, the Landau effective mass for protons 
and neutrons depends on the parameters Ci and Zj, where Ci is connected with the saturation 
properties of nuclear matter and Zj with the density dependence of the nuclear symmetry energy. 
Secondly, m* depends on the baryon density via the variable u, as well as on the isospin parameter 
/. Finally, there is also a direct dependence of m* on the regulator function g(k, Aj) via the 
parameters Aj. Hence, one can conclude that the values of the effective mass are reflected on the 
specific properties of the applied nuclear models. This enables us to carry out a systematic study 
of the properties of m* by applying the parametrization of the present model. 



3 Results and Discussion 

We have studied thermal and momentum-dependent effects on the properties of hot asymmetric 
nuclear matter including mean field potentials (proton, neutron and symmetry) as well as the 
effective mass, by applying a momentum-dependent effective interaction. The parametrization of 
the model can be found in Table 1 and Table 2 of Ref . pQ . 

In particular, we have studied two general cases. In the first one (called g{) the momentum 
regulator function g(fc,Aj) is given by equation (jSJ). In this case, the quantity Jl defined by 
Eq. ((7j), at T = 0, takes the form 

1 \ K f) \ W f W f 

where the upper signs refer to neutrons and the lower ones to protons. 

In the second case (called g-i) the momentum regulator function Aj) is given by equation 
(Q. Here, the finite range interactions are approximated as effective local interactions by retaining 
only the quadratic momentum dependence [U [14] . In this case, the energy density takes the form 
of Skyrme's effective interaction. At T = the quantity Jl is 





The parameters A, B, er, C\, Ci and B' (a small parameter introduced to maintain causality), 
are determined from constraints provided by the empirical properties of symmetric nuclear matter 
at the equilibrium density hq = 0.16 fm -3 . With the appropriate choice of the parameters, it is 
possible to parametrically vary the nuclear incompressibility K, so that the dependence on the 
stiffness of the equation of state may be explored. In the same way, by choosing the appropriate 
parameters x , x 3 , Z\ and Z 2 , it is possible to obtain different forms for the density dependence 
of the nuclear symmetry energy [fl [26] . 

Fig. 1(a) displays the behavior of the NSE as a function of the ratio u = n/n for the three 
different parametrizations of the function F(u) (relation [14]) . The function F^(u) leads to a stiffer 
nuclear symmetry energy dependence on the density, while the function F\ (u) leads to a softer one. 
It is worthwhile to point out that the above parametrization of the interacting part of the nuclear 
symmetry energy is extensively used for the study of neutron star properties [HES], as well as for 
the study of the collisions of neutron-rich heavy ions with intermediate energies [321 ES]- For a 
very recent review of the applications of the proposed momentum dependent effective interaction 
model and the specific parametrization of it, see Ref. [M] (and references therein). The aim of 
the above simple parametrization is to reproduce the nuclear symmetry energy originating from 
various theoretical calculations and/or experimental predictions and also to be able to cover the 
possible range of the nuclear symmetry energy dependence on the density. 

In Fig. 1(b) we plot the regulator momentum function g(k,A\) for the two cases applied in 
the present work. The most striking feature of the two curves is the fact that for low values of 
the momenta k, the two cases coincide, but for higher values of k they have completely different 
behavior. The motivation for applying these two different parameterizations for the momentum 
part of the interaction is to investigate in greater detail the momentum dependence effects on the 
mean field properties of hot asymmetric matter. 

The nuclear symmetry potential, pertaining to the isovector part of the nucleon mean field 
potential of a nucleon in nuclear matter, also depends on the momentum of the nucleon. In order 
to see the temperature and momentum effects on the nuclear symmetry potential, we first study 
the above effects on the nucleon single-particle potentials in hot asymmetric nuclear matter. 

In Fig. 2 we present the behavior of the proton and neutron single particle potentials for the 
case gi and three different values of the incompressibility K. The calculations are performed at 
the saturation density no- In each figure we display U p and U n for three different forms of the 
symmetry energy (named 1, 2 and 3 respectively, see Eq. [2]). It is noted that in case 3, for 
K = 180 MeV and for high values of k, U p and U n obtain comparable values. As we will see later, 
this affects the behavior of NSP. A feature of Fig. 2 is the effect of K on the values of the SPP 
for high values of k. Thus, for small values of k in the three cases, both for protons and neutrons, 
lead to similar values for the SPP, while for high values of k, there is an obvious splitting of the 
values of SPP, which shows a monotonic increase with K. 

Fig. 3 displays the trend of U p and U n , in case gi, for K = 240 and F 2 (u) = u as a function of 
the asymmetry parameter / = (n n — n p )/n. The meaning of the trend is very clear. The increase 
of the isospin asymmetry parameter / leads to a more attractive U p and also more repulsive U n . 
Similar behavior is found in the other cases as well. 

In Fig. 4 we plot U p>n for the case g 2 and for the three different values of K. The competitive 
behavior of U p and U n , both for low and high values, is clearly reflected in the behavior of U sym 
as indicated later in Fig. 7. 

The temperature effects on U p , U n are shown in Fig. 5. As expected, the increase of T leads to 
a corresponding increase of the values of U p and U n . The effects are more pronounced for T > 5 
MeV. The temperature effect on U p , U n has a significant influence on the temperature dependence 
of the nuclear symmetry energy as found in previous works [2TI [S] . 

Most of the theoretical models predict a decreasing symmetry potential with increasing nucleon 
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momentum, albeit at different rates, while a few nuclear effective interactions used in some models 
show opposite behavior. In the present work we try to clarify the above controversial results 
by applying a systematic study of U sym . In Fig. 6, we plot the nuclear symmetry potential 
U S ym = (U n — Up) 1 21 as a function of the nucleon kinetic energy E kin in case g±, for three values 
of K and in each case for three different choices of F(u). In the present work, we apply I = 0.4 
in the calculation of U sym and as pointed out also in US] , U sym is almost independent of the 
asymmetry parameter / but depends strongly on density and momentum. It is of interest to see 
that for the nine, in total, different cases, only in two of them U sym decreases with increasing 
nucleon momentum ( for K = 180 MeV and F^(u) = 2u 2 /(u + 1) and for K = 240 MeV and 
F(u) = u l/2 ). This trend is in agreement with experiment, which is well reproduced using the 
empirical relation ( 1221) . In the remaining seven cases, U sym increases with increasing nucleon 
momentum. In addition, from Fig. 6 it is concluded that the energy dependence of U sym is 
sensitive to the incompressibility K. 

It is worthwhile to notice that the nuclear symmetry potentials differs from the nuclear sym- 
metry energy as the latter involves the integration of the isospin-dependent mean-field potential of 
a nucleon over its momentum [8]. However, it is of interest to study the effect of the potential part 
of E sym {u) on the momentum and density dependence of U sym . Fig. 6 demonstrates the strong 
dependence of U sym on the function F(u). It seems that by applying the functions F{u) = u 1 ^ 2 
and F(u) = u we receive similar results only for low values of E kin (except for K = 120 MeV 
where a similar trend is taken also for higher values of E kin ). The case with F(u) = 2u 2 /(u + 1), 
which introduces a much stiffer density dependence E sym , leads to different results even for low 
values of momenta. To sum up, the density dependence of E sym is well reflected on the momentum 
dependence of U sym , so it is expected that models with different density dependence in the nuclear 
symmetry energy may predict different energy dependence of U sym . 

In Fig. 7 we also plot U sym as a function of E kin , for the case g 2 - In fact, the trends for K = 120 
MeV and K = 240 MeV coincide, and U sym exhibits a negative slope, while for K = 180 MeV the 
slope of U sym is positive. 

The key quantity to explain the trend of U sym as a function of E kin is expression (T2"4"|) . For 
this expression we conclude that the behavior of U sym (n, A;) as a function of k is defined from the 
values of the parameters C{ and Zj as well as from the function g(k,Ai). Both the parameters 
and Zi are related to the strength of the momentum dependence. In addition, Ci is fixed also to 
reproduce the properties of symmetric nuclear matter at the saturation point, while Zi is set to 
fix the density dependence of the symmetry energy. The regulator function g(k, Aj) mainly affects 
the trend of the slope of U sym . In particular, in case g 2 , after some algebra, the expression f[2"4l) 
can be written as 

U^ P (n, k) = U™{n, a, Z t ) - V 1 E km u £ (<7 '~ 8Z<) , (31) 

i=l,2 A i 

where V 1 is a constant. Eq. (jHTl) establishes a linear relation between U sym and E kin , in case g 2 , 
as indicated in Fig. 9. 

In case g\ for low values of k, a similar relation as fl3Tl) holds. Accordingly, for high values of 
k we find 

U^ P (n, k)^^uY: {d - SZ t ) A 2 , (32) 

I^kin 1=1,2 

where T> 2 is a constant. 

From the above analysis it is clear that the behavior of U sym depends on the specific values of 
the parameters Cj, Zi and Aj and also on combinations of them. In particular, the values of the 
parameters define the sign of the slope (positive or negative) and the regulator function g(k,Ai) 
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defines the trend of the slope. Thus, the role of the regulator function g(k, A$) is of importance, 
concerning the momentum-dependent behavior of the proton and neutron single particle potential, 
as well as the symmetry potential. 

In view of the above comment it may be interesting to explore other choices for the regulator 
function as those used in Ref. |65j. Such work is in progress. 

It is worth noting that according to Fig. 7, U sym is strongly related to the density dependence 
of the nuclear symmetry energy and the total parametrization of the interaction energy density 
part, i.e. for various values of the incompressibility K. 

The three cases where U sym is a decreasing function of nucleon momentum are shown in Fig. 8 
and compared with the phenomenological expression (122]) . 

In Fig. 9 we plot U sym as a function of E^ in for the nuclear symmetry energy parametrization 
Fi(u) = w 1 / 2 and F 2 (u) = u and for three different values of u. In both cases, U sym , exhibits strong 
density dependent behavior. In particular, there is a monotonic increase of U sym as function of 
the nucleon density. 

Special effort has been devoted to the study of the thermal effects on U sym . So far, that 
problem has received little theoretical attention [8]. Temperature influences the first term of the 
right-hand side of Eq. (1251) i.e. U^^f^n, I). The thermal effects on U sym are displayed in Fig. 10. 
In the first case (Fig. 10[a]) an increase of T leads to decreasing values of U sym , while the inverse 
behavior is seen in the second case (Fig. 10 [b]). However, in both cases thermal effects do play a 
role, especially for T > 5. Temperature does not change the trend, but only slightly affects the 
values of U sym depending on Ej~ in . 

Another quantity of interest, which can be easily calculated, is the isoscalar potential. Several 
decades ago, it was already pointed out that the quantity {U n + U p )/2, which is obviously the 
single-nucleon potential in absence of asymmetry, should be a reasonable approximation to the 
isoscalar part of the optical potential. The momentum dependence of (U n + U p )/2 is important for 
extracting information about the symmetric matter equation of state [4T] • Now, in order to check 
the validity of the model parameters this is customary, and a more stringent test to compare the 
isoscalar potentials, as calculated in the present work, with them of the variational many-body 
(VMB) predictions by Wiringa [531157]. 

In Fig. 11 we plot the isoscalar potential at the four values of momenta k (k = 1, 2, 3, 4 fm _1 ) 
for four different cases. The results of the present work compared with those (named UV14+UVII) 
predicted by Wiringa [51]. In Fig. 11(a) indicated that for the specific case, the predictions of the 
present work are in very good agreement with the VMB predictions up to about k = 3 fm _1 . In a 
second case, presented in Fig. 11(b), there is agreement only for low values of the density and for 
momenta k. There is obvious disagreement for the predictions presented in Figs. 11(c) and 11(d). 
From the above analysis it is concluded that by suitable choice of the model parameters, we 
are able to reproduce and/or be consistent with the predictions of other microscopic many-body 
calculations or experimental constraints. This is an advantage of the present model that is flexible 
enough to reproduce predictions of many other theoretical models. 

The Landau effective mass splitting is shown in Fig. 12, where the effective masses of proton 
and neutron are displayed as functions of the asymmetry parameter / for all cases corresponding 
to g\. For all cases, m* is a decreasing function of /, while m* is an increasing function of /. 
The only exception is the case for K = 240 MeV with Fs(u) = 2u 2 /(u + 1), where m* is a 
decreasing function of /. In every case the isospin mass splitting (m* — m*) is positive and is 
directly dependent on the incompressibility K and the parametrization function F(u). 

In the present work, we also address the problem which arises when the effective mass splitting 
directly affects the energy dependence of the symmetry potential [32]. We found that in the 
framework of the proposed model and for the regulator function gx, all the cases lead to a positive 
splitting but the slope of the U sym may be positive or negative. So, it is concluded that the trend 
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of U sym , by applying the present parametrization, does not directly connect to the effective mass 
splitting. Nonetheless, we found that there is a correlation between U sym and m* — m* originated 
from the density dependence of E sym . More precisely, the cases with F(u) = u 1 ^ 2 and F(u) = u 
which lead to much closer values for U sym also lead to a similar value of m* — m* as is displayed 
in Figs. 12 and 13. 

To illustrate further the dependence of the effective mass on the asymmetry parameter / and 
to find the quantitative characteristic of this dependence, the values of m*/m for various values 
of / have been derived with the least-squares fit method and found to take the general form 

m*(I) 

~ c + cj + c 2 J 2 . (33) 

m 

The parameters Cj depend on the model parameters Ci, Zi,Ai and they are different for each 
case. Specifically, for low values of / it is obvious that linear relations between m*/m and / is hold 
(see also Fig. 12). The splitting of effective masses, according to Eq. (13"3"|) . is well approximated 
by 

Am* {I) = m* n (I) - m* p {I) ~ dj + d 2 I 2 , (34) 

where the parameters di also depend on the model parameters Ci, Zi, A*. For low values of J, the 
linear relation Am* (I) m d\I is hold. 

The effective mass as a function of u is indicated in Fig. 13. In all cases the effective mass 
splitting (m* — m*) is positive. 

In Fig. 14(a) the effective mass, for the case g 2 , is displayed. The splitting (m* — m*) is positive 
for the case with K = 120 MeV (and also for K = 240 MeV) and negative for K = 180 MeV. 
Thus, one can conclude by also comparing with the case gi, that the regulator function g(k, Aj) 
can have a dramatic effect on the effective mass splitting. The above feature is displayed also 
in Fig. 14(b), where the effective mass is plotted as a function of u. It is worth pointing out 
that in case g 2 , an almost linear relation between m*/m and / is hold even for high values of I 
(in comparison with the case gi). So, the relation (|33|) with c 2 = 0, clearly describes the above 
dependence. 

It should be emphasized that according to Figs. 12 and 13, the density dependence of the 
symmetry energy influences the values and also the slope of the effective mass (both for proton 
and neutron), but does not affect the sign of the mass splitting, which always remains positive 
(case gi). Nevertheless, in case g 2 , for different values of K the sign of the mass splitting may 
be negative or positive. In the case with K — 180 the positive slope of U sym (see Fig. 7) is 
connected with the negative splitting (see Fig. 14 (a)). This result is in accordance with the finding 
of Ref. [371 E2]. 

From the above analysis one can conclude that, in general, and by applying the proposed, in 
the present momentum effective interaction model, the slope (positive or negative ) of the U sym as 
a function of Eki n is not connected directly with the sign of the effective mass splitting. Certainly, 
there is a connection between the above quantities in the sense that both, according to Eqs. ([2"7]) . 
(I28p and (j24p . depend in a similar way on the parameters Ci, Zi Aj. 



4 Summary 

In this study we have applied a momentum dependent effective interaction in order to investi- 
gate the single-particle properties of hot asymmetric nuclear matter. More specifically, we have 
examined the single-particle potentials of protons and neutrons, the asymmetry potential and the 
isospin mass splitting for various cases. The effects of the specific parametrization of the inter- 
action part of the energy on the single-particle properties are studied and analyzed. It has been 
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concluded that the behavior of the symmetry potential depends strongly on the parametrization of 
the interaction part of the energy density as well as on the momentum dependence of the regulator 
function. The effects of the parametrization are less pronounced on the isospin mass splitting. 
The effect of an increase of the temperature is just to shift higher the values of the proton and 
neutron single particle potential. The symmetry potential U sym can be an increasing or decreasing 
function of the nucleon kinetic energy, depending on the parametrization of the momentum de- 
pendent effective interaction model. In the first case U sym increases with the temperature, while 
in the second one the inverse behavior is observed. 
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Figure 1: (a) The nuclear symmetry energy as a function of u = u/uq for three different 
parametrizations of the function F(u) given by Eq. (TH1) . (b) The regulator momentum function 
g(k,A) for the two cases, given by Eqs. ([8]) and ([9]), applied in the present work. For more details 
see text. 




Figure 2: The proton and neutron single particle potentials (the upper curves corresponds to U n , 
while the lower one to U p ) for the case gi and for three different vales of the incompressibility K 
(at the saturation density uq). In each figure we display U p and U n for three different forms of the 
nuclear symmetry energy (named 1,2 and 3 respectively). 
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Figure 3: U p and U n , in case g\ and for K = 240 MeV and i^u) = u versus k and for various 
values of the asymmetry parameter I = [n n — n p )/n. 




Figure 4: U p (left) and U n (right) for the case #2 and for three different values of K (at the 
saturation density n ). 
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Figure 5: U p and U n , in case g\ and for K = 240 and F 2 {u) = u as a function of temperature T 
(in MeV) (at the saturation density n ). 




E kjn (MeV) E kn (MeV) E kin (MeV) 



Figure 6: The nuclear symmetry potential U sym = {U n — U P )/2I versus the nucleon kinetic energy 
Ekm m case gi, for three values of K and in each case for three different choices of F(u) (at the 
saturation density n ). 
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Figure 7: U sym = {U n — U P )/2I versus Ekin in case g%, for three values of K (the cases K = 120 
MeV and K = 240 MeV are coincide) (at the saturation density n ). 




Figure 8: U sym = (U n — U P )/2I versus E kin for the three cases (at the saturation density n ), 
which decreases with increasing of nucleon momentum, in comparison with U sym constrained by 
the experimental data. For more details see text. 
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Figure 9: U sym = (U n — U P )/2I versus Ef-m in case g%, for K = 240 MeV for three different values 
of u (a) for Fi(u) = u 1 ^ 2 and (b) for F 2 (u) = u. 
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Figure 10: U sym = {U n — U p )/2I versus E kin in case gi, for K = 240 MeV for various values of T 
(in MeV) (a) for F^u) = u 1 ' 2 and (b) for F 2 {u) = u. 
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Figure 11: The isoscalar potential (U n + U p ) /2, as a function of baryon density, at the four values of 
the momenta k, for four different cases, in comparison with the variational many-body calculations 
(UV14+UVII). 




Figure 12: The effective mass of proton and neutron as a function of the asymmetry parameter / 
for all the cases corresponding to g\ (at the saturation density n ). 



21 




1 2 3 4 5 










(c) 




9i 




K=240 MeV 




1=0.4 




- P 2 




- "—Pa 





























1 2 3 4 5 



$ure 13: The effective mass of proton and neutron versus u for all the cases corresponding 



0.9 



0.8 



0.7 



0.6 



I I 1 1 



(a) 
9 2 

F 2 (u)=u 

— ■ — p (K=120) 
— • — n (K=120) 
■ ■ ■ ■ p (K=180) 
• ■ •■ ■ n (K=180) 



_i i i i i i i_ 



0.0 0.2 0.4 



p. 6 0.8 1.0 



1.0 



0.8 



0.6 



0.4 



0.2 



1 1 1 1 



(b) 




B 2 




1=0.4 






p (K=120) 




n (K=120) 




p (K=180) 




n (K=180) 




_l 1 I 1 I 1 l_ 



2 u 3 



$ure 14: (a) The effective masses of proton and neutron versus I for two cases corresponding 
(b) The effective masses of proton and neutron versus u for two cases corresponding to gi. 
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